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Abstract 

The Daya Bay and RENO reactor neutrino experiments have revealed that the small- 
est neutrino mixing angle is in fact relatively large, i.e., # 13 ~ 9°. Motivated by this 
exciting progress, we perform a systematic study of the neutrino mass matrix M v with 
one or two texture zeros, in the assumption that neutrinos are Dirac particles. Among 
fifteen possible patterns with two texture zeros, only three turn out to be favored by 
current neutrino oscillation data at the 3<r level. Although all the six patterns with one 
texture zero are compatible with the experimental data at the 3a level, the parameter 
space of each pattern is strictly constrained. Phenomenological implications of M v on 
the leptonic CP violation and neutrino mass spectrum are explored, and the stability of 
texture zeros against the radiative corrections is also discussed. 
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1 Introduction 



The solar, atmospheric, accelerator and reactor neutrino experiments have provided us with 
compelling evidence that neutrinos are massive particles and they can transform from one 
flavor to another [1]. The lepton flavor mixing can be described by a 3 x 3 unitary matrix 
U, which is usually parameterized through three flavor mixing angles (^12^23^13) an d one 
CP- violating phase 5. To be explicit, we adopt the following parametrization 

/ C 12 C 13 S 12 C 13 S 13 \ 

U = 



C 12 S 23 S 13 S 12 C 23 e %& S 12 S 23 S 13 + C 12 C 23 e ^ S 23 C 13 

\ — CinCnoSio + SioSnoe * — Si 9 Cqo Si q — Ci n Sqq Q % CnqCiq/ 



(1) 



-12 u 23°13 ~ °12°23' z - c, 12^23' : '13 u 12 23° ^23^13- 

where s^ = sin#^ and = cos^ (for ij = 12, 23, 13) have been defined. If neutrinos 
are Majorana particles, two additional CP-violating phases (p, a) have to be introduced to 
fully describe the flavor mixing. Thanks to the elegant neutrino oscillation experiments, two 
neutrino mass-squared differences (5m 2 , Am 2 ) and two flavor mixing angles (#12,^23) have 
been measured with a reasonably good precision. More recently, the Daya Bay [2] and RENO 
[3] collaborations have clearly observed the disappearance of v e from nuclear reactors, and 
revealed that the smallest mixing angle is relatively large, i.e., 9 13 ~ 9°. This is really a great 
news to the long-baseline neutrino oscillation experiments, which aim to pin down the sign of 
Am 2 and the magnitude of the CP-violating phase 5. In spite of the great progress made in 
neutrino physics, our understanding of neutrino properties is far from complete. For instance, 
the absolute scale of neutrino masses is not yet determined and whether neutrinos are Dirac 
or Majorana particles remains an open question. 

Since a convincing flavor theory is lacking, the approach of texture zeros has been suggested 
to study the flavor problem for a long time [4]. The texture zeros of a fermion mass matrix 
dynamically mean that the corresponding matrix elements are sufficiently suppressed in com- 
parison with their neighboring counterparts [5], and they can help us to establish some simple 
and testable relations between flavor mixing angles and fermion mass ratios [6]. In fact, a great 
number of works have been devoted to confronting the zero textures of neutrino mass matrix 
with the neutrino oscillation data [7, 8, 9, 10, 11, 12, 13, 14]. However, most of them have 
assumed neutrinos to be Majorana particles, because various seesaw mechanisms for neutrino 
mass generation lead to light Majorana neutrinos. 

If neutrinos are Dirac particles, they can acquire masses exactly in the same way as quarks 
and charged leptons do in the standard model. In this scenario, it seems quite difficult to ex- 
plain why the neutrino Yukawa couplings are twelve orders of magnitude smaller than the top- 
quark Yukawa coupling. Although this mass hierarchy problem has never been well understood 
even for charged fermions, it has been shown that the highly-suppressed Yukawa couplings for 
Dirac neutrinos can naturally be achieved in the models with extra spacial dimensions [15] or 
through radiative mechanisms [16]. Conservatively speaking, the most important motivation 
for considering Dirac neutrinos is the fact that no experiments have already excluded such a 
possibility. In the present work, we simply assume neutrinos to be Dirac particles and perform 
a systematic study of the neutrino mass matrix M v with one or two texture zeros. 
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Without loss of generality, we can take the mass matrix M v for Dirac neutrinos to be Her- 
mitian by redefining the right-handed neutrino fields. As M v is Hermitian, three independent 
off-diagonal matrix elements are in general complex, while three independent diagonal ones are 
real. If n of them are taken to be vanishing (i.e., M v has n independent texture zeros), then 
we shall arrive at 



n\ (6 — n)\ 



(2) 



different textures. There are totally fifteen two-zero textures of M u , which can be classified 
into six categories: 





t ° 





A 







X 


A 




\A* 


A* 


X 




( x 


A 





Bi = 


A* 





A 




v o 


A* 


X 




f x 





A 


B 3 : 








A 




VA* 


A* 


X 



( 

A* 

V o 



A 
x 

A* 



0\ 

A 
x/ 





( X 







B 2 : 





X 


') 




VA* 


A* 


/ 




( X 


A 




B 4 : 


A* 


X 


: ) 




V o 


A* 


/ 



(3) 



(4) 



C : 




(5) 



Di 



I X 

A* 

\A* 



A A\ 

x / 



( X 

A* 

\A* 



A A\ 
x 
/ 



(6) 



E i: 



and 



Fx: 



( ° 

A* 

\A* 



A 


A* 



A\ 
A 
x / 




E, 



2 • 



( 

A* 
\A* 



A A\ 
x A 
A* / 




E, 



/ A A\ 

A* x 
\A* x/ 




(7) 



(8) 



in which each "x" stands for a nonzero and real matrix element, while each "A" for a nonzero 
and complex one. Note that this classification is similar to that for the two-zero textures of 
Majorana neutrino mass matrix [7, 8, 9, 10], which are symmetric and complex rather than 
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Hermitian. Although the one-zero textures are in general less predictive than the two- zero ones, 
we shall consider them for completeness. Assuming one of six independent matrix elements to 
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where the notations are the same as those in Eqs. (3)-(8). It is straightforward to observe 
that three two-zero patterns F 123 in Eq. (8) can be excluded, because they lead to just one 
nonzero flavor mixing angle. 

We aim to confront the remaining twelve two-zero textures and six one-zero textures with 
the latest global-fit results of current neutrino oscillation data done by Fogli et al. [17] 1 . Only 
three two-zero textures, i.e., A l5 A 2 , and C, are found to be compatible with current experi- 
mental data at the 3a level, so are all the six one-zero patterns, i.e., P i (for i — 1, 2, • • • , 6). In 
particular, most physical consequences of those viable patterns with one or two texture zeros 
have been explored in an analytical way, and the stability of texture zeros against radiative 
corrections is also discussed. We establish the relationship between the location of texture 
zeros and CP conservation, and demonstrate that the two- zero patterns A 12 and the one-zero 
patterns P 4)5i6 lead to CP conservation in the lepton sector. 

The remaining parts of this paper are organized as follows. In section 2, we give some 
general remarks on the texture zeros for Dirac neutrinos. We show that it is possible to fully 
determine the neutrino mass spectrum and the CP-violating phase 5 for all two-zero textures, 
as well as for the one-zero textures P 



4,5,6- 



The relationship between the location of texture 
zeros and CP violation is pointed out. The stability of texture zeros for Dirac neutrino mass 
matrix against the renormalization group running is considered. Section 3 is devoted to the 
analytical and numerical analyses of the two-zero textures of M v , while section 4 to the one-zero 
textures. Finally we summarize our conclusions in section 5. 



2 General Remarks 
2.1 Important relations 

In the flavor basis where the charged-lepton mass matrix M l is diagonal, the Dirac neutrino 
mass matrix M v can be reconstructed in terms of three neutrino masses (m 1 ,m 2 ,m 3 ) and the 

lr rhe global-fit analysis of neutrino oscillation experiments has also been done in Refs. [18, 19], however, 
all the results are completely consistent with each other at the 3cr level. Therefore, our discussions will not be 
affected when the different global-fit results are used. 
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flavor mixing matrix U. Namely, 

M„ = U\ A, \l\ (11) 




where A x = 77 • m 1 , A 2 = x ' m 2 an d ^3 = m 3 with 77, x = ±1. Note that the three eigenvalues of 
a general 3x3 Hermitian matrix are real, but not necessarily positive, so we have chosen the 
signs of the first two eigenvalues relative to the third one as rj and x- The parametrization of 
U through three flavor mixing angles (d 12) 6 23 ,6 13 ) an d one CP-violating phase 5 has already 
been given in Eq. (1). 

Now we explore the general consequences of the texture zeros of M v . If one element of 
M v is vanishing, i.e., (M u ) aj3 = 0, then we can obtain the corresponding constraint among the 
flavor mixing parameters 

v ■ m.u^u;, + x ■ m 2 u a2 u; 2 + m 3 u a3 u; 3 = . (12) 

Note that Eq. (12) implies one constraint condition in the case of a — f3 (e.g., the one-zero 
textures Pi 2,3)) but two constraint conditions in the case of a 7^ j3 (e.g., the one-zero textures 
P 4 5 6 ). In the former case, one can derive from Eq. (12) that either f] = x = — 1 or i] ■ x = — 1 
must hold and 

where £ = m 1 /m 3 and ( = m 2 /m 3 have been defined. In the latter case, we can get two 
constraint conditions by requiring both real and imaginary parts of the left-hand side of Eq. 
(12) to be zero. More explicitly, we have 

(itf - xC) • Im [Kg] = , (14) 

and 



Re 



-^23 



i-xC 



-"■13 



(15) 



where K?f = U^U^U^U^j. As will be shown later, the relations in Eqs. (13), (14) and (15) 
are very useful in the determination of flavor mixing parameters when we discuss the one- zero 
textures in section 4. 

If two independent elements of M u are vanishing, i.e., (M v ) ab = (M v ) a/3 = with ab 7^ a(3 
as shown in Eqs. (3)- (8), we can obtain 

^~m 3 V ' U a2 U* b2 U al U* pi -U al U* bl U a2 Ufc' 

; "^3 " x ' u a2 u* b2 u al u* pi - u al u* bl u a2 u* p2 ■ 1 } 

Since both ^ and ( are by definition real and non-negative, the imaginary parts of the quantities 
on the right-hand side of Eq. (16) have to disappear. This requirement may lead us to the 
determination of the CP-violating phase S, as we shall show below. 
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2.2 Texture zeros and CP violation 



Now that 6 13 has been measured to be relatively large, the CP-violating effects are promising 
to be discovered in the long-baseline neutrino oscillation experiments, if the CP-violating phase 
5 turns out to be not extremely small or very close to ir. As is well known, the CP violation 
is characterized by the Jarlskog invariant J, which is defined as [20] 



Im 



K 



7 k 



(17) 



where e aj3 ^ and e ijk denote the Levi-Civita symbol, and K°f has been defined below Eq. (15). 
It is now straightforward to discuss the CP violation for the zero textures. Some comments 
are in order: 



CP-violating one-zero textures - In the case of (M u 
either 77^ = xC 01 im 



K. 



a 8 
23 



= with a 7^ /3, Eq. (14) implies 
±J~ = 0. Note that the solar neutrino experiments have 
established m 2 > m 1 , or equivalently ( > £, so we are left with J = 0. Therefore, CP 
violation is only possible for the patterns Pi i2 ,3, and we have 5 = or n for the other 
one-zero patterns. 

• CP-violating two-zero textures - Note that Eqs. (13), (14) and (15) apply as well to the 
two-zero textures. Hence CP violation is only possible for the textures with both a = b 
and a = p. In other words, the patterns C, E l5 E 2 can lead to leptonic CP violation, 
while 5 = or n holds for all the other two- zero textures. 

In order to demonstrate the above observation of CP-violating two-zero textures, we can 
directly calculate the imaginary part of the first identity in Eq. (16). It turns out that 
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K ah 
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+ (ab <-> a(3) = , (18) 



where u ab a (3" stands for the foregoing terms with the exchange of a -h- a and b -h- j3 both 
in the superscripts and in the subscripts. Since the patterns F 12S have already been excluded, 
we need to consider only two different cases: (1) a = b and a = (3; (2) a = b and a 7^ j3. 
In the first case, it is easy to see that the quantities on the right-hand side of Eq. (16) are 
automatically real, so there is no constraint on the CP-violating phase. In the second case, 
Eq. (18) reduces to Im K^f = ±J = 0, which is consistent with our previous observation. 
Note that the result for the second identity in Eq. (16) can be obtained by exchanging the 
subscripts "1" and "2" in Eq. (18), however, it doesn't give any new constraints. Hence we can 
conclude that only one off-diagonal texture zero in the Hermitian Dirac neutrino mass matrix 
is enough to ensure CP conservation in the lepton sector. 



2.3 Parameter counting 

Since we have assumed massive neutrinos to be Dirac particles, there are seven physical pa- 
rameters: three neutrino masses (m 1 ,m 2 , m 3 ), three flavor mixing angles (9 12 , 9 23 , 9 13 ), and one 
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CP-violating phase 5. The number of constraint relations caused by the texture zeros depends 
on the location of zeros, so we count the parameters for four distinct situations: 



1. One off-diagonal zero - This category contains the patterns P 4i5i 6- As shown in Eqs. 
(14) and (15), there are two independent constraint relations. Therefore, with the help 
of current experimental measurements of three flavor mixing angles (6* 12 , # 23 , 9 13 ) and two 
neutrino mass-squared differences, defined as [17] 

5m 2 = m\ — m\ , Am 2 = ml — - (m\ + m£) , (19) 

we can fully determine the neutrino mass spectrum and the CP-violating phase. As 
shown in section 2.2, these patterns predict CP conservation and thus 5 = or n. Using 
Eq. (15) and the following relation 

R _ a™ 2 _ 2(c 2 -a (2(]) 

u ~\Am 2 \ |2-(C 2 + e)l' 
we can fix £ and (. Thus three neutrino masses are given by 

m 3 = ^2 _ ^2 ' m 2 = m sC , m 1 = m 3 £ . (21) 
The detailed discussions of P 4i5i6 will be given in section 4. 

2. One diagonal zero - This category contains the patterns P x 2 3- The texture zero leads 
to one constraint condition in Eq. (13). Given the five experimental observables, one 
parameter is left free for these patterns. If the CP-violating phase is fixed, we can pin 
down the neutrino mass spectrum by using Eqs. (13) and (20). On the other hand, if 
the lightest neutrino mass is assumed, one can determine the CP-violating phase. 

3. One diagonal zero and one off-diagonal zero - This category contains the patterns A 12 , 
B lj2i3j4 , D 12 , and E 3 . These two texture zeros impose three constraint conditions, 
so both three neutrino masses and the CP-violating phase can be calculated by using 
experimental observables. Since these patterns are CP conserving, we have 5 = or ir. 
Therefore, £ and ( can be calculated from Eq. (16), and the neutrino mass spectrum is 
then given by Eq. (21). Note that £ and £ have to satisfy Eq. (20), so there will be one 
testable correlative relation among the flavor mixing angles and neutrino masses. 

4. Two diagonal zeros - These category contains the patterns C and E 12 . The texture 
zeros induce two constraint relations, which together with five experimental observables 
leads to the full determination of neutrino mass spectrum and the CP-violating phase. 
This can be done as follows. First, note that £ and ( are functions of the CP- violating 
phase 5 as shown in Eq. (16). Then it is possible to determine or constrain 5 from Eq. 
(20). Once 5 is fixed, we can obtain (£, () from Eq. (16), and thus neutrino masses from 
Eq. (21). All the two-zero textures will be discussed in great detail in section 3. 
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Now we summarize the latest global-fit results of three flavor mixing angles (9 12 , 9 23 , 9 l3 ) 
and two neutrino mass-squared differences (5m 2 , Am 2 ), which will be taken as experimental 
observables to determine the other flavor parameters. For the normal mass hierarchy with 
Am 2 > 0, it has been found at the 3a level [17] 

0.259 < sin 2 9 12 < 0.359 or 30.6° < 9 12 < 36.8° , 
0.331 < sin 2 9 23 < 0.637 or 35.1° < 9 23 < 53.0° , 
0.017 < sin 2 13 < 0.031 or 7.5° < 9 13 < 10.1° ; (22) 

and 

6.99 x 10~ 5 eV 2 < 5m 2 < 8.18 x 10" 5 eV 2 , 

2.19 x 10~ 3 eV 2 < +Am 2 < 2.62 x 10~ 3 eV 2 . (23) 

For the inverted mass hierarchy with Am 2 < 0, the global analysis yields [17] 

0.259 < sin 2 9 12 < 0.359 or 30.6° < 9 12 < 36.8° , 
0.335 < sin 2 9 23 < 0.663 or 35.4° < 9 23 < 54.5° , 
0.017 < sin 2 9 13 < 0.032 or 7.5° < 9 13 < 10.3° ; (24) 

and 

6.99 x 10" 5 eV 2 < 5m 2 < 8.18 x 10" 5 eV 2 , 

2.17 x 10" 3 eV 2 < -Am 2 < 2.61 x 10~ 3 eV 2 , (25) 

at the 3a level. The best-fit values of three mixing angles are 9 l2 = 33.6°, 9 23 = 38.4°, and 
9 13 = 8.9°, while those of neutrino mass-squared differences are 5m 2 = 7.54 x 10" 5 eV 2 and 
Am 2 = 2.43 x 10~ 3 eV 2 . It is interesting to note that the best-fit value of the CP-violating 
phase is 5 ~ n, which happens to be consistent with the prediction of the Dirac mass matrix 
with one off-diagonal texture zero. However, there is no constraint on 5 at the 3a level. 



2.4 Stability of texture zeros 



The stability of texture zeros for Majorana neutrinos has already been discussed in the litera- 
ture [21, 12], by using the renormalizat ion-group equations (RGEs) [22]. It has been demon- 
strated that the texture zeros of Majorana neutrino mass matrix are stable against the one-loop 
quantum corrections [12]. In this subsection, we shall examine whether the stability of texture 
zeros for Dirac neutrinos is maintained [14]. 

To accommodate Dirac neutrino masses, one can extend the standard model with three 
right-handed neutrino singlets, and simply require the lepton number conservation to forbid 
the Majorana neutrino mass term. At the one-loop level, the RGEs for Dirac neutrinos and 
charged leptons can be written as [23] 



16tt 



16tt 



dt 
dt 



a u + CI (YX) + C l u (Ytf) 



Y, 



(26) 



where t = \n(/j,/M z ) with fi being an arbitrary renormalization scale, and M z the Z-boson 
mass. Here Y ul = \f2M vl /v denotes respectively the neutrino and charged- lept on Yukawa 
coupling matrices, and v ~ 246 GeV is the vacuum expectation value of the Higgs field. 

In the standard model, the relevant coefficients are = C\ = +3/2, C l v = C\ = —3/2, 
a u = — 9gf/20 — 9g|/4 + T and a { = — 9gf/4 — 9gf/4 + T, where g 1 and g 2 are gauge couplings, 
and T = Tr [3 (l^Y?) + 3 (Y d Y d f ) + (Y u Yj) + (y z Y/)] with Y q (for q = u, d) being the up- 
type and down-type quark Yukawa coupling matrices. The term Y u Yj in Eq. (26) can be 
safely neglected, because neutrino masses are much smaller than charged-lepton masses. In 
the flavor basis where the charged-lepton Yukawa matrix is diagonal Y { = Diag{y e , y M , y T }, one 
can observe from the second identity in Eq. (26) that the charged-lepton Yukawa coupling 
matrix remains diagonal as the energy scale evolves. After solving the RGE for Y u , one can 
find that the neutrino mass matrix at the weak scale M z is related to that at a high-energy 
scale A in the following way 



M V {M Z ) = J 



(h 


Vo 








M„(A) , 



(27) 



where the RGE evolution function I represents the overall contribution from gauge and quark 
Yukawa couplings, while I a (for a — e, //, r) stand for the contributions from charged-lepton 
Yukawa couplings. Because of m e m T , we have I e < 1^ < I T and they will modify 

the structure of M v . In contrast, J 7^ 1 just affects the absolute scale of neutrino masses. 
However, the texture zeros of M v are stable against the one-loop RGE corrections. Taking the 
pattern A-,^ for example, we have 





a* 




(28) 



at A, and thus 








al e 
c*L dL 



(29) 



at M z . Note that the neutrino matrix M v 1 {M Z ) is no longer exactly Hermitian. Nevertheless, 
it can be shown that I a ~ 1 (for a = e, /i, r) hold as an excellent approximation in the 
standard model. As a consequence, the previous important relations derived from texture 
zeros are formally valid both at A and M z , at least in the lowest-order approximation. 



3 Two-zero Textures 

3.1 Analytical approximations 

First of all, we consider the two-zero textures and explore their phenomenological implications 
in an analytical way. It has been demonstrated that a permutation symmetry between the 
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viable two- zero patterns A x and A 2 , B x and B 2 , or B 3 and B 4 for Majorana neutrinos [12]. 
The existence of such a symmetry originates from the fact that the location of texture zeros 
in each pair is related by an exchange between the last two rows and columns of M u . For 
the same reason, it is only necessary to study the patterns A l5 B x , B 3 , D 1; E x and E 3 of 
Dirac neutrino mass matrix, and the results for the other patterns can be obtained with the 
replacements 9 23 — > n/2 — 9 23 and 5 — > 5 — ir. Now we discuss these patterns by making some 
reasonable approximations. 



Pattern A x with (M u ) ee = {M v ) = 0. As we have discussed, there is no CP violation 
and thus 5 = or ir. In the case of 5 — 0, we obtain from Eq. (16) that 



£ = +v 



£l3 / s 12 s 23 _ 

c 13 \ c 12 c 23 / 



C = -X-^l^ + s*) ■ (30) 



2 

C 13 \ S 12 C 23 



Obviously, r] = +1 and x = — 1 should be taken to ensure that £ and ( are non- negative. 
Since s^ 3 1 still holds, Eq. (30) approximates to 

£ ps sin 9 13 tan 9 23 tan , 

( ~ sin tan # 23 cot . (31) 

Given 0.59 < tan^ 12 < 0.75, 0.70 < tan^ 23 < 1.3, and 0.13 < sin^ 13 < 0.18 at the 
3a level, Eq. (31) leads to £ < ( < 1. Hence only the normal neutrino mass hierarchy 
is allowed, i.e., Am 2 > 0. Taking the best-fit values of three neutrino mixing angles 
(i.e., 9 12 = 33.6°, # 23 = 38.4°, and 9 13 = 8.9°) and those of two neutrino mass-squared 
differences (i.e., 5m 2 = 7.54 x 10~ 5 eV 2 and Am 2 = 2.43 x 10" 3 eV 2 ), one can figure out 
the neutrino mass spectrum 



m 3 



VAm 2 = 4.9 x 10" 2 eV , 
m 2 ~ m 3 sin 13 tan0 23 cot0 12 = 9.0 x 10~ 3 eV , 

m 3 sin 9 13 tan 9 23 tan 9 12 = 4.0 x 10~ 3 eV . (32) 



m 1 

As shown in Eq. (20), there exists an interesting correlation between three mixing angles 
and the neutrino mass ratios 

= Sm^ _ 4 tan 2 ^ 23 sin 2 9 13 
v Am 2 ~ sin2# 12 tan2# 12 ' 1 ' 

Taking the best-fit values of three mixing angles, we can get R v ps 0.027 from Eq. (33). 
On the other hand, the best-fit values of two neutrino mass-squared differences yield 
R v ~ 0.031. Therefore, Pattern A-,^ is well consistent with current oscillation data, and 
will be soon tested in the future neutrino oscillation experiments. In the case of 5 = ir, 
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we can obtain 

c 13 \ c 12 c 23 / 
c 13 \ 6 12 c 23 / 

where 77 = —1 and x = +1 are implied. To the leading order of s 13 , the phenomenological 
implications in the case of 5 = n are the same as those in the case of 5 = 0, so we shall 
not discuss this case further. More precise measurements of neutrino mixing angles are 
needed to distinguish between these two cases. 

Pattern A 2 with (M v ) ee = {M v ) eT = 0. All the analytical results of Pattern A 2 can be 
obtained from those of Pattern A x with the replacements 9 23 — > n/2 — 9 23 and 5 — > 5 — n. 
Therefore, in the case of 5 = 0, we have 

c 13 \ c 12 s 23 / 

c = +x .^(^l_ s \ . (35) 

c 13 \ 6 12 6 23 / 

After taking 77 = — 1 and x — +1 an d neglecting the terms of 0(s\ 3 ), one can obtain 

£ ps sin # 13 cot tan # 12 , 

£ sin 9 13 cot # 23 cot 9 l2 • (36) 

Given 0.59 < tan# 12 < 0.75, 0.70 < tan# 23 < 1.3, and 0.13 < sin# 13 < 0.18 at the 
3cr level, we can verify that £ < £ < 1, implying that only the normal neutrino mass 
hierarchy is allowed. The neutrino mass spectrum turns out to be 



m 3 vArn? , 

m 2 ~ m 3 sin # 13 cot # 2 3 tan # 12 , 

mj ~ m 3 sin 6* 13 cot 9 23 cot # 12 . (37) 

As in the case of Pattern A x , there is a constraint relation among three mixing angles 
and neutrino mass ratios. To the leading order, we get 

5m 2 ^ 4 cot 2 # 23 sin 2 9 13 
v = Ari 2 ~ sin 29 12 tan 29 12 ' ( ' 

Although Eq. (38) is not fulfilled by the best-fit values of (9 12 , 9 23 , 9 13 ) and (5m 2 , Am 2 ), 
Pattern A 2 is indeed compatible with current oscillation data at the 3a level, as will 
be demonstrated by the numerical analysis in subsection 3.2. In a similar way, one can 
discuss the case of 5 = n, for which the analytical results can be obtained from Eqs. (30) 
and (31) by replacing 9 23 with n/2 — 9 23 . 
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Pattern B x with (M„) w = (Mj eT = 0. With the help of Eq. (16), we obtain 

£ _ _ ~ S 12 C 12 g 23 C 13 ~F ( C 12 C 23 + g 12 g 23) C 23 S 13 + 2s 12 C 12 S 23 C 23 S 2 3 
S 12 C 12 S 23 C 23 ^ ( S 12 — C 12) C 23 S 13 + S 12 C 12 S 23 S 13(1 + C 23) 

£ _ ^, _ ~ S 12 C 12 S 23 C 13 ^ ( C 12 C 23 + S 12 S 23) C 23 S 13 + ^ S 12 C 12 g 23 C 23 S 13 ^gg^ 
S 12 C 12 S 23 C 23 ^ ( S 12 — C 12) C 23 S 13 + S 12 C 12 S 23 S 13(1 + C 23) 

where the upper and lower signs refer to the cases of 5 = and 5 = ir, respectively. In 
the leading order approximation, one can get i] = \ = ~~ 1 an d £ ~ C ~ tan 2 # 23 . In the 
next-to-leading order approximation, we find 

4sin# 13 

£ - C ~ ± — • (40) 

sin 26 12 sin 2# 23 v 1 

Since 5m 2 > or equivalently ( > £, only <5 = 7r is allowed. In this case, the constraint 
relation turns out to be 

R >*W^t Ml - (41) 

Note that R v oc sin# 13 in Eq. (41) may be one order of magnitude larger than R v oc 
sin 2 9 13 in Eq. (33) or Eq. (38). Taking the values of (9 12 , 9 23 , 9 13 ) m the 3a ranges, one 
verify that R v > 0.8, which is obviously in conflict with the experimental observation 
R v < 0.038. Therefore, we conclude that Pattern B x has already been excluded by cur- 
rent neutrino oscillation data. Due to the permutation symmetry between Pattern B x 
and Pattern B 2 , Eq. (41) is also applicable to the latter, implying that Pattern B 2 
has been experimentally ruled out as well. 



Pattern B 3 with (MJ W = {M v ) e(l = 0. With the help of Eq. (16), we obtain 

£ = -v 



1 fifi \ u/efj, 

S 23 S 12 S 23 ~F C 12 C 23 S 13 



c 23 s 12 c 23 =t C 12 S 23 S 13 
£ — . £23 . C 12 g 23 ^ S 12 C 23 S 13 ^) 

C 23 C 12 C 23 "F S 12 S 23 S 13 

where the upper and lower signs refer to the cases of S = and 5 = ir, respectively. In the 
leading order approximation, rj = \ — ~ 1 an d £ ~ C ~ tan 2 9 23 hold, as for Pattern B x . 
However, in the next-to-leading order approximation, one gets 

4 tan 2 L sin 9, , , 

£-C~T ■ ofl 23 ofl 13 , (43) 
sin 29 12 sin 29 23 

and 

2 sin 2 

^ ~ • n 13 tan 2 # 23 1 tan 29 23 1 . (44) 
sin ZD 

It is evident from Eq. (43) that only 5 = is allowed, because of ( > £. By taking the 
values of (9 12 , 9 23 , 9 l3 ) in the 3a ranges, we can obtain R v > 0.4, which is far outside of the 
3a range of R v . Hence Pattern B 3 is not compatible with current neutrino oscillation 
data, nor is Pattern B 4 due to the permutation symmetry. 
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Pattern C with (MJ^ = (MJ TT = 0. With the help of Eq. (16), we obtain 

c 12 c 13 c 12( c 23 — S 23) ~~ C ^ s l2 s 23 c 23 s l?, c 8 



i = v 



'13 25x20x28230230,5(1 + sf 3 ) — (cf 2 — Sl2)( c 23 — s 23) s 13 



( — X ' Sl2ClS • s 12( s 23 c 23) 2c 12 g23 c 23 g 13 c 3 

2-Sl2 c 12 s 23 c 23 c <5(l + S 13 ) — (cf 2 — Sl2)( c 23 — s 23) s 13 



'13 



where we have defined c s = cos 5. Note that the leptonic CP violation is allowed in this 
case. Generally speaking, it is inappropriate to expand Eq. (45) in terms of s 13 and ignore 
the higher-order terms, because c s and c 2 3 — S23 could be vanishingly small according to 
current neutrino oscillation data. But the analytical approximations in some interesting 
limits deserve further discussions: 

1. 5 — tt/2 and #23 7^ vr/4. Insetting 5 = tt/2 into Eq. (45), one arrives at £ = 
—T] cot 2 9 13 cos 2 9 12 / cos29 12 and ( = +x cot 2 9 13 sin 2 9 12 / cos 26> 12 , which are inde- 
pendent of #23- Since cos 2 9 12 > sin 2 9 12 holds, we are led to rj = —1, x — +1> 
and £ > C- This observation indicates that the maximal CP-violating phase is not 
allowed if 9 23 7^ it/ 4. 

2. 5 7^ it/2 and 9 23 = 7r/4. Assuming #23 = vr/4 in Eq. (45), we can get r\ = x = — 1 
and £ = C = cos2 ^13/(1 + sin 2 6> 13 ) < 1, which only depends on 9 13 . Hence the 
maximal mixing 9 23 — 7r /4 is not favored if 5 7^ ir/2, and a tiny deviation of 9 23 
from 7r/4 is necessary to break the degeneracy between m l and m 2 . 

3. 5 = it/2 and 9 23 = it/A. In this case, one can verify that {U^ 2 = \U Ti \ 2 (for 
% = 1, 2, 3) hold, so Eq. (16) is not applicable. The equality (M„) w = (M V ) TT implies 
that only one constraint relation is obtained by requiring (M v ) w = 0. Assuming 
5 = tt/2 and 6*23 = tt/A, we find 

^cj 2 sj 3 + s\ 2 ) + C(s 2 2 s 2 3 + c\ 2 ) = c 2 3 (46) 

in the case of rj = x = — 1- If the higher-order terms of 0(s\ 3 ) are neglected, then 
one obtains £ = (1 — ( ■ c 2 2 )/s 2 2 , implying that neither ( > £ > 1 nor 1 > £ > £ is 
allowed. Consequently, both £ and £ should be close to one, and the deviations from 
one are of 0(s\ 3 ). Based on this observation, we define £ = 1 — A£ and ( = 1 — A(, 
and then solve Eqs. (20) and (46) for A£ and A(. It is straightforward to get 



A£ 



(2 + it>Jsin 2 fl 13 
1- R u cos29 12 /2 ' 



AC = (2 ~^ )sin2 ^ (47) 
S l-i^cos20 12 /2 " 1 ; 

Hence the neutrino masses are nearly degenerate, while both 5 ~ 7r/2 and # 2 3 ~ 7r /4 
are allowed and expected to be valid simultaneously. 
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4. 5 ^ tt/2 and 9 23 ^ ir/4. Furthermore, we assume c 5 ^> s 13 and then obtain 

cot L cot 29 



/ 12 z^ 23 

sin # 13 cos 5 



. / tan 6* 12 cot 2# 23 \ 

C ~ -X • 1 + — f • 48 

\ sin 6> 13 cos o J 

If cot 2# 23 cos 5 > holds, one can set x — ~ 1 an d then obtain ( > 1, indicating 
the inverted mass hierarchy m 2 > m 1 > m 3 . In this case, £ is required to be larger 
than one but smaller than £, and this can be achieved by setting rj — +1 and 
cot 12 cot 2^23 > 2 sin 13 cos S > 2 cot 2# 12 cot 2# 23 . It is straightforward to verify 
that cot 2# 23 cos 5 < contradicts with the requirement ( > £. Therefore, only the 
inverted mass hierarchy is allowed. 

From the above discussions, one can observe that there exists a small parameter space 
around 5 = n/2 and #23 = vr/4, where the normal mass hierarchy is allowed and neutrino 
masses are in fact nearly degenerate. In the main part of parameter space, only the 
inverted mass hierarchy is consistent with current oscillation data. For the general case, 
the numerical analysis will be done in subsection 3.2. 

Pattern D x with (MJ = (M„) = 0. With the help of Eq. (16), we obtain 



£ w cot# 12 tan# 23 , 

( ps tan 9 12 tan 9 23 , (49) 

where 77 = —1 and x = +1 are taken, and the higher-order terms 0(s\ 3 ) have been 
safely neglected. Together with the oscillation data, Eq. (49) indicates ( < £, which is 
in contradiction with the experimental observation m 2 > m 1 . Therefore, this pattern is 
not viable, nor is Pattern D 2 due to the permutation symmetry. 

Pattern E 1 with (Mj ee = (Mj w = 0. With the help of Eq. (16), we obtain 

£ ps sin 2 9 12 tan 2 9 23 j c °s 2# 12 , 

C Ps cos 2 6> 12 tan 2 6> 23 / cos 2# 12 , (50) 

in the leading-order approximation, where we have set rj = +1 and x = — 1. Since the 
neutrino mass ratios have been determined as in Eq. (50), one can get 

R „ 4tan 4 fl 23 cos2fl 12 

"~ |tan 4 ^ 23 (l + cos 2 2^ 12 ) -4cos 2 2# 12 | ' 1 ; 

implying R u > 2cos2# 12 or R u > tan 4 # 23 / cos2# 12 . Taking the allowed values of 9 23 and 
9 l2 in the 3a ranges, we can derive R v > 0.56 or R u > 0.50, respectively. Therefore, 
Pattern E x is not consistent with current neutrino oscillation data, nor is Pattern E 2 
according to the permutation symmetry. 
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• Pattern E 3 with (Mj ee = (MJ MT = 0. With the help of Eq. (16), we obtain 

f « sin 2 fl 12 /cos2# 12 , 

C « cos 2 # 12 /cos2# 12 , (52) 

in the leading-order approximation, where we have chosen rj = — 1 and x = +1. The 
ratio of neutrino mass-squared differences turns out to be 

which is obviously disfavored by current experimental data. 

In summary, we conclude that only three two-zero patterns (i.e., A 1; A 2 and C) are 
consistent with current neutrino oscillation data. The detailed numerical analysis of these 
three patterns will be performed in the following subsection. 



3.2 Numerical Analysis 

We have numerically confirmed that only three patterns A x , A 2 and C are viable, as obtained 
from the above analytical approximations. Our strategy for numerical calculations is as follows: 

1. For each of the twelve patterns of M v in Eqs. (3)-(7) we generate a set of random numbers 
of (9 12 , 9 23 , 9 13 , 5m 2 ) lying in their 3a ranges, which are already given in Eqs. (22)-(25). 
For the CP-conserving patterns, we consider both cases of 5 = and 5 = re. For the 
CP- violating patterns, 5 is allowed to vary in the range [0, 2ir). 

2. With the above generated numbers, we can calculate the other physical parameters of 
M v . First of all, the neutrino mass ratios £ and ( can be determined through three 
mixing angles and the CP-violating phase. To be consistent with the experimental data, 
the following two conditions should be satisfied: (a) m 2 > m 1 or equivalently ( 2 > £ 2 ; (b) 
since only the neutrino mass hierarchies m 2 > m 1 > m 3 and m 3 > m 2 > m 1 are allowed, 
we further require (( 2 — 1)(£ 2 — 1) > 0. With the calculated neutrino mass ratios £ and 
C, one can figure out the absolute neutrino masses m i (for i = 1,2,3) by inserting the 
randomly generated 5m 2 into Eq. (21). Finally we compute Am 2 via neutrino masses 
m iy the patterns can be regarded as viable if the computed Am 2 is consistent with the 
experimental values in Eqs. (23) and (25). 

3. From all the points satisfying the above consistency conditions, we can figure out three 
neutrino mass eigenvalues (m 1 ,m 2 ,m 3 ) via Eq. (21). For the CP-violating patterns, it is 
possible to calculate the Jarlskog invariant J = Si2 c i2 s 23 c 23 s i3 c i3 s <5- To present numerical 
results, we show the allowed regions of three neutrino mass eigenvalues (m 1 ,m 2 ,m 3 ), 
those of three mixing angles (9 12 , 9 23 ) 9 13 ), and those of (J, 5) if CP violation is allowed. 

Our numerical results are shown in Figs. 1-3. Some comments and discussions are in order: 
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• Pattern A x - Our numerical results show that the patterns A x and A 2 survive current 
experimental tests, which confirms our analytical calculations. The allowed ranges of 
neutrino mass eigenvalues and mixing angles for Pattern A t and Pattern A 2 are de- 
picted in Fig. 1 and Fig. 2, respectively. Since the patterns A x and A 2 are similar to 
each other due to the permutation symmetry, here we focus on the former. From Fig. 
1, some interesting observations should be noted: (a) Only the normal mass hierarchy 
m 1 < m 2 < m 3 is allowed, (b) The deviation of 9 23 from 7r/4 is expected in both cases of 
5 = and 5 = n. However, 6> 23 < 45° is favored for 5 = 0, while 9 23 > 45° for 5 = n. (c) 
The parameter space in the case of 5 = receives more stringent constraints than that 
in the case of 5 — tt. Therefore, more precise measurements of neutrino mixing angles 
are required to distinguish between the cases of 5 = and 5 = ir, as well as between 
Pattern A 1 and Pattern A 2 . 

• Pattern C Unlike the previous two patterns, the leptonic CP violation is allowed 
for this pattern. Our numerical results are given in Fig. 3, where only the inverted 
mass hierarchy has been considered. From the analytical analysis in section 3.1, the 
normal mass hierarchy is indeed allowed, if 6* 23 ~ 7r/4 and 5 ~ ir/2. However, the main 
parameter space points to the inverted neutrino mass hierarchy. Although the mixing 
angles for this pattern are not strictly constrained, the discovery of CP violation in the 
future long-baseline neutrino oscillation experiments will definitely single out Pattern C 
among all the two- zero textures as the true Dirac neutrino mass matrix. 

Finally it is worth pointing out that we have performed a numerical analysis of the two- zero 
textures of M v by the oscillation data at the 2a level. We find that all these three patterns are 
compatible with current experimental data at this level, although the corresponding parameter 
space is somewhat smaller. 

4 One-zero Textures 

Now we proceed to consider the one-zero textures of Dirac neutrino mass matrix M v . Since 
only one zero element is assumed, these textures are not as predictive as the two-zero ones. On 
the other hand, it will be no doubt that all the six patterns in Eqs. (9) and (10) are compatible 
with current experimental data, which has also been confirmed by numerical calculations. For 
simplicity, we focus on the analytical analysis by making some reasonable approximations. 

• Pattern V x with (Mj ee = 0. With the help of Eq. (12), we get 

V ■ &l 2 c\ z + X ■ (s 2 12 c n + 4 3 = • ( 54 ) 

Since only one constraint relation exists, it is by no means possible to pin down both the 
neutrino mass spectrum and the CP-violating phase. From Eq. (54), one can observe 
that 5 is entirely unconstrained. In this case, however, neutrino mass eigenvalues can 
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be determined, as we shall show below. To find the solutions to Eq. (54), one has to 
consider three different possibilities: (1) rj = \ = ~ lj (2) V = — 1 and x = +lj (3) 
f] = +1 and x — ~ 1- Now we examine whether all these possibilities are allowed by 
current experimental data. 

1. If i] = x = — 1 is assumed, one can observe from Eq. (54) that both £ and £ should 
be of 0{s\ 3 ), and furthermore obtain 

£ = tan 2 13 /sin 2 £ 12 -£cot 2 # 12 , (55) 

from which £ < tan 2 9 13 can be derived by requiring £ > £. Since both £ and £ are 
quite small in this case, R v ps £ 2 — £ 2 holds as an excellent approximation. To be 
explicit, one can find 



tan 2 9 13 cot 2 9 12 — y tan 4 9 13 + R v cos 2# 12 
^ ~ cot 2 9 12 - 1 



tan 2 13 tan 2 12 — J tan 4 13 + R v cos 2# 12 
C = tan 2 fl 12 -l ' (56) 

which together with Eq. (21) leads to the neutrino mass spectrum. In addition, the 
condition tan 2 9 l3 > y^R^sin 2 9 l2 has to be fulfilled to guarantee a non- negative £. 
Given 0.259 < sin 2 9 12 < 0.359, 0.017 < tan 2 9 13 < 0.032, and 0.027 < R v < 0.037 
at the 3(7 level for the normal neutrino mass hierarchy, it is straightforward to verify 
that such a condition cannot be satisfied. Hence there is no solution in this case. 

2. If 7] — +1 and x — ~ 1 are taken, then Eq. (54) can be written as 

£ = tan 2 fl 13 /sin : % + £cot 2 # 12 ■ (57) 

Different from the previous case, both £ and £ need not to be as small as tan 2 # 13 . 
Requiring £ > £ leads to cot 2 9 12 > 1, which is favored by current oscillation data. 
Inserting Eq. (57) into Eq. (20), we can figure out the neutrino mass ratios 



^ tan 2 9 13 cot 2 9 12 - ^tan 4 9 13 + R u cos 2# 12 



1 - cot 2 9 12 



tan 2 6-, o tan 2 9, 9 + A /tan 4 0, , + R cos 29 



- 13 lml 1/12 -T V 13 ' v 12 , s 

C " l-tan 2 12 ' (58) 

where we have assumed £ 2 < £ 2 1, and thus R v ~ £ 2 — £ 2 . For the best-fit values 
12 = 33.6°, 13 = 8.9°, and R v = 0.031, the neutrino mass ratios are given as £ = 
0.045 and £ = 0.22, which justifies the assumption of £ 2 < £ 2 1. With the help 



of Eq. (21), the neutrino mass eigenvalues are found to be m 3 ~ y5m 2 / (£ 2 — £ 2 ) 
0.04 eV, m 2 = (m 3 ps 8.8 x 10 -3 eV and m l = £m 3 m 1.8 x 10 -3 eV. 
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3. If i] = — 1 and x — +1 are taken, then Eq. (54) turns out to be 

( = £ cot 2 9 12 - tan 2 9 13 / sin 2 12 . (59) 

The requirement ( > £ gives rise to a lower bound £ > tan 2 # 13 / cos26 ) 12 - In a 
similar way to the previous case, we assume £ 2 < ( <C 1 and thus i2 v ~ ( 2 — £ 2 . 
Consequently, we arrive at 

tan 2 13 + tan 2 9 12 \J tan 4 6 ) 13 + R v cos 2# 12 
^ = 1 - tan 2 12 ' 

tan 2 9 13 tan 2 12 + Vtan 4 9 13 + R u cos2^ 12 
C = l-tan 2 12 • (60) 

For the best-fit values 9 12 = 33.6°, 9 13 = 8.9°, and R v = 0.031, the neutrino 
mass ratios are given as £ = 0.13 and ( = 0.22, which justifies the assumption of 
£ 2 < C 2 <^ 1- The neutrino masses turn out to be m 3 = 5m 2 / (( 2 — £ 2 ) ps 0.049 eV, 
m 2 = £m 3 0.01 eV, and m l = £m 3 ps 6.4 x 10~ 3 eV. 

So we conclude that Pattern P 1 is consistent with current neutrino oscillation data, 
and only the normal mass hierarchy is allowed. Additionally, the CP-violating phase 5 
is arbitrary. 

Pattern P 2 with (MJ^ = 0. With the help of Eq. (12), we get 

r ?'^( s 12 c 23 + 2j7"t^ C 13 +C 12 S 23 S 13 ) + X'C( C 12 C 23~ ^"^^8 C 13 + S 12 S 23 S 13) + S 23 C 13 = ) (61) 

where J = Si2 c i2 s 23 c 23 s i3 c i3 s <5 1S the Jarlskog invariant in the standard parametrization, 
and t s = tan 5 has been defined. For simplicity, we assume 77 = x = — 1 and solve Eq. 
(61) for £ and ( in the leading-order approximation. In this case, it is straightforward to 
observe £ ~ ( ~ tan 2 9 23 . In order to figure out the deviations of £ and ( from tan 2 9 23 , 
we define £ = tan 2 9 23 — A£ and ( = tan 2 9 23 — A(, and insert them back into Eq. (61). 
Then one can see that both A£ and A( should be of 0(s\ 3 ). Combining Eq. (61) with 
Eq. (20), we obtain 

_ (1 - tan 4 9 23 )R U + 2 tan 4 9 23 sin 2 fl 13 /(cos 2 6 12 cos 2 9 23 ) 
? ~ 2tan 2 23 /cos 2 12 -tan 2 23 (l -tan 2 ^)^ ' 1 ' 

and A( = sin 2 9 13 tan 2 9 23 /(cos 2 9 l2 cos 2 9 23 ) — A£tan 2 ^ 12 , where we have assumed the 
normal mass hierarchy. Taking the best-fit values of (0 12 , 9 23 ) 9 13 ) and (5m 2 , Am 2 ) for 
example, one can get £ = 0.59 and ( = 0.61, and the neutrino masses are m 3 = 6.4 x 
10~ 2 eV, m 2 = 3.9 x 10~ 2 eV and m 1 = 3.8 x 10 -2 eV. Whether the neutrino mass 
hierarchy is normal or inverted depends on the mixing angle 9 23 . Similarly, one can 
also analyze the case with r] ■ x — ~ 1- It is worthwhile to point out that 5 is almost 
irrelevant to the determination of neutrino masses, because its contribution to Eq. (61) 
is suppressed by sin0 13 . Therefore, 5 is totally arbitrary, as for Pattern P x . 
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Because of the permutation symmetry, the analytical results for Pattern P 3 can be 
obtained by the replacements 9 23 — > n/2 — 9 23 and 5 — > n — 5. Thus we shall not discuss 
further about this case. 

Pattern P 4 with (Mj eM = 0. With the help of Eq. (12), we get 

~ 1 1 ' £ C 12( S 12 C 23 ^ C 12 S 23 S 13) + X ' C S 12( C 12 C 23 "F S 12 S 23 S 13) + S 13 S 23 = ? (63) 

where the upper and lower sign stands for 5 = and 5 = n, respectively. At the leading 
order, Eq. (63) approximates to (r?£ — xC) = 2 sin 9 13 tan 9 23 j sin2# 12 , implying x = —\ 
and 77 = =bl. Take x — — 1 and rj = +1 for example. Since = 2 sin 9 13 tan 9 23 / sin 29 12 , 
we see that R v ( 2 — £ 2 holds as an excellent approximation. Consequently, one gets 



4 tan 2 6> 23 sin 2 6> 13 - it!,, sin 2 29 12 
4 sin 26 l 12 tan 9 23 sin # 13 



= 4 tan 2 9 23 sin 2 fl 13 + R v sin 2 2fl 12 
4 sin 29 12 tan 6* 23 sin 13 

implying the normal neutrino mass hierarchy. Taking the best-fit values of (9 12 , 9 23 , 9 13 ) 
and (Sm 2 , Am 2 ), we have £ = 0.07 and ( = 0.19, and thus the neutrino masses m 3 = 
4.9 x 10" 2 eV, m 2 = 9.3 x 10~ 3 eV and m 1 = 3.4 x 10~ 3 eV. In the case x = -1 and 
77 = —1, it is straightforward to verify that only the inverted mass hierarchy is allowed. 

Because of the permutation symmetry, the analytical results for Pattern P 5 can be 
obtained by the replacements 9 23 — > n/2 — 9 23 and 5 — > n — 5. 

Pattern P 6 : (M v ) = 0, we have 
r]i[s\ 2 =F 2 cot 29 23 c l2 s l2 s l3 - c 2 12 s 2 13 ] + xC[c\ 2 ± 2 cot 29 23 c l2 s l2 s l3 - s 2 l2 s 2 l3 ] = c\ 3 , (65) 

where the upper and lower sign stands for 5 = and 5 = ir, respectively. In the leading- 
order approximation, we can obtain rj^s\ 2 + xC c i2 = c i3- Note that neutrino oscillation 
data indicate c 2 2 > s\ 2 and ( > £, so only two possibilities need to be considered: (1) 
77 = x — +lj and (2) 77 = — 1 and x — +1- Now we examine whether these two 
possibilities are compatible with current neutrino oscillation data: 

1. If 77 = % = +1 is assumed, then £ £ cos 2 # 13 holds. In order to figure out the 
deviations of £ and £ from cos 2 13 , we define £ = cos 2 13 — A£ and £ = cos 2 9 l3 — A(, 
and insert them back into Eq. (65). Then we obtain A( = — tan 2 12 A£. Combing 
Eq. (65) with Eq. (20), we obtain A£ ps R v sin 2 9 13 cos 2 9 12 . Taking the best-fit 
values of neutrino mixing parameters, we finally get m 3 = 0.38 eV. Considering 
m 1 ps m 2 ps m 3 , we can see this might violate the cosmological limit J2 m i = 
m x + m 2 + m 3 < 0.61 eV from the cosmic microwave background and the large-scale 
structure observations [24]. 
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2. If i] = — 1 and x = +1 are taken, we can get £s 2 2 + 1 ~ C c i2- Current oscilla- 
tion data indicate s\ 2 < c 2 2 , thus £ > 1 is required, implying the inverted neu- 
trino mass hierarchy m 3 < m l < m 2 . Furthermore, combining £ ps £cot 2 # 12 — 
l/sin 2 fl 1 2 > 1 with 2(C 2 - £ 2 )/(C 2 + £ 2 ), we can obtain C = l/(cos 2 # 12 - 

sin 2 12 ^J(2 + R u )/{2 — R v )). Taking the best-fit values of neutrino mixing pa- 
rameters, one can obtain ( = 2.62 and £ = 2.60, and the neutrino masses are 
m 3 = 2.7 x 10~ 2 eV, m 2 = 7.1 x 10~ 2 eV and m 1 = 7.0 x 10~ 2 eV. 

Therefore, we conclude that Pattern P 6 is completely consistent with current neutrino 
oscillation data, and only the inverted mass hierarchy is allowed. However, the absolute 
neutrino masses might exceed the cosmological bound in the case of a nearly-degenerate 
neutrino mass spectrum. 

5 Summary 

In light of the recent measurements of 9 13 and the latest global-fit analysis of neutrino oscillation 
experiments, we have performed a systematic study of the Dirac neutrino mass matrix M v 
with two independent texture zeros or one texture zero. It turns out that three two-zero 
patterns (i.e., A 12 and C) and all six one- zero patterns(i.e., Pi 5 2,3,4,5,6) can survive current 
experimental tests at the 3a level. In fact, all of them are found to be compatible with the 
oscillation data even at the 2a level. Analytical analyses have been done for both two-zero 
and one-zero textures in order to understand why they are favored or disfavored by current 
experimental data. Moreover, the allowed parameter space for three viable two- zero textures 
has been obtained through a detailed numerical analysis. The following is a brief summary of 
our conclusions: 

• In the basis where the charged-lepton mass matrix is diagonal, neutrino masses, flavor 
mixing angles, and CP-violating phase are determined by the Dirac neutrino mass matrix 
M u , which can be further made Hermitian by redefining the right-handed singlet neutrino 
fields. We demonstrate that one vanishing off-diagonal element in the Hermitian neutrino 
mass matrix is enough to guarantee CP conservation in the lepton sector. 

• Among fifteen two-zero textures of M v , only three (i.e., A l5 A 2 and C) turn out to be 
consistent with current neutrino oscillation data at the 3a level. We have explained in 
detail why the other patterns are disfavored. For Pattern A 1 and Pattern A 2 , only the 
normal neutrino mass hierarchy is allowed, and there is no CP violation. For Pattern C, 
both normal and inverted mass hierarchies are in principle allowed. However, the normal 
mass hierarchy is possible only if 5 n/2 and 6> 23 7r/4. The precise measurements of 
neutrino mixing angles and leptonic CP violation are needed to distinguish among these 
currently viable patterns. 

• All the six one-zero textures are compatible with current neutrino oscillation data. For 
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the patterns P 4 5 6 , even only one texture zero is assumed, it is possible to fully deter- 
mine neutrino mass eigenvalues and the CP-violating phase in terms of the observed 
three mixing angles (9 l2 , #23> ^13) an d two neutrino mass-squared differences (5m 2 , Am 2 ). 
For the patterns P 1 2 3 , even if one free parameter remains, the full determination of neu- 
trino mass spectrum is possible, at least in the first-order approximation. Although the 
parameter space for all these patterns has received some constraints, it seems impossible 
to exclude them experimentally in the near future. 

It is worthwhile to stress that whether neutrinos are Dirac or Majorana particles remains 
an open question. At present, the only feasible way to demonstrate that neutrinos are Ma- 
jorana particles is to observe the neutrinoless double-beta decays [25]. If neutrinos are Dirac 
particles, such a lepton-number-violating process is forbidden. Therefore, it really makes sense 
to consider Dirac neutrinos, and to study the flavor mixing and CP violation in this scenario. 
The ongoing and upcoming neutrino oscillation experiments are expected to precisely measure 
the neutrino mixing parameters, in particular the smallest mixing angle 9 13 , the deviation of 
6> 23 from 7r/4 and the Dirac CP-violating phase 5. The sensitivity of future cosmological obser- 
vations to the sum of neutrino masses J2 m i wm probably reach ~ 0.05 eV in the near future. 
We therefore expect that some patterns of the two-zero textures of the neutrino mass matrix 
M v for Dirac neutrinos might be excluded or only marginally allowed by tomorrow's data, and 
those surviving the experimental tests should shed light on the underlying flavor structure of 
massive neutrinos. 
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Figure 1: Pattern A x of M v : The allowed regions of neutrino masses (m l , m 2 , m 3 ) and neu- 
trino mixing angles (sin 2 9 12 , sin 2 9 23 , sin 2 9 13 ), where the left column is for 5 = while the right 
column for 5 = n. The global-fit data of two neutrino mass-squared differences (5m 2 , Am 2 ) 
and three neutrino mixing angles (0 12 , 23 , 8 13 ) at the 3a level have been input. 
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Figure 2: Pattern A 2 of M u : The allowed regions of neutrino masses (m l , m 2 , m 3 ) and neu- 
trino mixing angles (sin 2 9 12 , sin 2 9 23 , sin 2 9 13 ), where the left column is for 5 = while the right 
column for 5 = n. The global-fit data of two neutrino mass-squared differences (5m 2 , Am 2 ) 
and three neutrino mixing angles (0 12 , 23 , 8 13 ) at the 3a level have been input. 
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Figure 3: Pattern C of M u : The allowed regions of neutrino masses (m 1 ,m 2 ,m 3 ), neutrino 
mixing angles (sin 2 9 12 , sin 2 9 23 , sin 2 9 13 ), and the Jarlskog invariant versus the CP-violating 
phase (J, 5). The global-fit data of two neutrino mass-squared differences (5m 2 , Am 2 ) and three 
neutrino mixing angles (^12^235^13) a ^ t ne 3a level have been input, while the CP-violating 
phase 5 is allowed to freely vary in [0, 2tt). 
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